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1. Introduction 

With the recent observation of Bose-Einstein conden- 
sate (BEC) in both systems of trapped rubidium [1] and 
sodium [2] atoms, and the reported evidence of BEC in 
a system of lithium atoms [3], it has became urgent to 
study the diffraction patterns in the near-resonant imag- 
ing of cold atomic samples. This problem is especially 
important since it concerns the only straight forward 
in situ detecting scheme for the tiny atomic clouds un- 
dergoing evaporation cooling [4] . Optical imaging pro- 
vides thus a possible diagnostics for the formation of 
BEC. 

In the experiment at JILA by Anderson et al. [1], a 
weakly interacting BEC (according to the criteria set by 



the BEC research in homogeneous systems [5]) has 
been observed. Even though such a condensate is closer 
to an ideal noninteracting BEC than any other previous 
system which exhibits bosonic degeneracy, its kinetic 
properties differ significantly from the ideal gas case. It 
is now well established that for the particular state of 
^^Rb used at JILA the scattering length is positive [6] 
(i.e., atom-atom interactions are predominantly repul- 
sive), while that of ^Li used by Bradley et al. is negative 
[7] (i.e., atom-atom interactions are attractive). The in- 
homogeneity introduced by a trap potential allows one 
to optimize the compression in phase space [8], but it 
strongly affects the boundary conditions for the system. 
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This is one of the reasons why BEC in a trap is possible 
for atoms with a negative scattering length [9], whereas 
it is forbidden in a homogeneous system [10]. Indeed, in 
a trap the relation between stability of the BEC and the 
two-body interaction scattering length becomes more 
complex even at zero temperature [9,11,12]. 

To probe the presence of BEC and to study its unusual 
properties, we have to consider its interactions with the 
external world. Scattering of photons from a near reso- 
nant laser provides one of the most natural choices for 
such a measurement scheme [13]. 

The aim of this paper is to examine the problem of 
imaging BEC within a cold atomic sample by coherent 
scattering of a weak, near-resonant, continuous-wave 
laser field. Two separate studies are required for the full 
description of near-resonant imaging. First, a reliable 
theory of the weak field scattering has to be developed, 
and the near field pattern has to be calculated. Second, 
realistic calculation of the image formation process is 
needed. Such a calculation must necessarily take into 
account the presence of apertures, lenses etc. In this 
paper we combine the recently developed theory of co- 
herent light scattering [13], with the theory of optical 
imaging. 

The paper is arranged as follows: in Sec. 2 we briefly 
review the theories of weak light scattering off a BEC. 
We concentrate on the two approaches developed in Ref 
[13]: the on-shell approximation, and Glauber's gener- 
alized diffraction theory. Using these methods, we con- 
struct the near field distribution of the scattered light, 
and express the scattering off the condensate in terms of 
an optical transmission process with a complex trans- 
mission function including both phase shift (refraction) 
and amplitude attenuation (absorption). We compare 
the numerical results for this transmission function ob- 
tained using the two above mentioned approaches. In 
Sec. 3 we describe a simple imaging system consisting 
of a single ideal converging lens and an aperture in front 
of it. The impulse response function is derived in the 
paraxial approximation limit, and the expression for the 
field distribution in the image plane is given. This is the 
main result of this work. In Sec. 4 we analyse numeri- 
cally the properties of the image using the parameters 
that correspond closely to those used in the JILA and 
Rice experiments [1,3]. Finally, we conclude in Sec. 5. 



and dense condensates should lead to formation of 
polaritons , a gap in the excitation spectrum, and deflec- 
tion, or even back reflection of the laser light from the 
atomic sample [14, 15]. Small condensates, on the other 
hand, should lead to a collective response of atoms, and 
a collectively broadened Lorentzian spectrum [16]. For 
not too large, and not too dense condensates (i.e., for 
situations corresponding to experimental conditions of 
Refs. [1, 3]), the theory developed by us in Refs. [13, 
17], seems to be the most suitable. This theory does not 
lead to the formation of a gap in the spectrum, since the 
photon energy is larger than the appropriately defined 
optical potential. Photons are thus scattered primarily in 
the forward direction, but nevertheless are affected 
strongly by the medium. The scattering and absorption 
cross sections have nonLorentzian shapes [17], with an 
overall width determined by collective excitations (as in 
Ref. [16]), 



7eff ^3Nyl2{Ka)\ 



(1) 



where A^ is the number of atoms, y is a single atom 
spontaneous emission width, ^l = 27r/A is the laser 
wavevector, X is the laser wavelength, and a is the con- 
densate size. The spectrum gradually narrows towards 
the resonance, and exhibits a cusp at the exact resonance 
with a width controlled by single atom dephasing pro- 
cesses and losses, such as spontaneous emission to 
uncondensed states, quantum diffusion etc. 

Using the expressions derived in Ref. [13], we can 
construct the near field distribution. In the on-shell 
approximation [18], our starting point is a scattering 
equation valid in the weak field limit. For spherically 
symmetric density distributions, the asymptotic solution 
for the averaged photon operator au^j, has the form 



+ B{k, iJL)d{ck - oyO^Q-'^^\ 



(2) 



where the first term describes the incident plane wave, 
while the second describes the process of elastic scatter- 
ing with the scattering amplitude 



2. Coherent Light Scattering 



B{k, M) = S f2Bi(h)Pi(cos0)(2U 1)7477, (3) 



The literature on coherent light scattering from cold 
atomic clouds is already quite extensive (see Ref. [13], 
and references therein). Qualitatively, various effects are 
predicted, depending on the size and density of the 
atomic cloud, in particular, of the condensate. Large 



where Pi (x) are the Legendre polynomials. We use here 
the coordinate system in which k = ki (sin cos <^, sin ^ 
sin <^, cos 0), Ul = Ul (0, 0, 1), and the polarization vector 
s^i^ = (1,0, 0), so that z defines the optical axis. The 
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coefficients ^/fe) can be well approximated by the 
expression 



Biih) 



jL 



[r+ y[ - i (wl - (Oo)] 



(4) 



i 



where y[ = lliry r^drn(r)j](kLrX with n(r) being 



the density distribution of cold atoms in a trap and ji 
being the spherical Bessel function of order /. F is an 
effective single atom dephasing rate which accounts for 
spontaneous emission to the uncondensed states, quan- 
tum diffusion of the excited wave packet, etc. Evidently, 
in typical situation F is larger than both the quantum 
diffusion rate and y[13]. (Ol and C(>o are respectively the 
probe laser frequency and the atomic resonance 
frequency. 

The reconstruction of the field distribution at any 
spatial point can be made according to 



Eosir) 



2 dkg(k) Sf,akf,d^ 



^o(p)e''^^tAos(p,zX 



(5) 



In the above expression Soip) is the probe field ampli- 
tude, whereas Q{k) is the dipole coupling constant, 
which is a slowly varying function of k and is related to 
y by y= (87r^^L/3c)|p(^L)2|l We assume that the po- 
larization of the field Eosir) is approximately equal to 
the polarization of the probe field, so vector notation can 
be suppressed. 

The scaled near field amplitude is given by 



ijJosip.z) =1+2 ^^BdkO 

1=0 ^ 



/; 



X sineded^^'^^'''^-^^Jo(ki^psme)Pi(cose), (6) 



where Jq is the Bessel function of order 0. In the gener- 
alized diffraction theory (GDT) [19], the solution ob- 
tained after the propagation through the cold atom sys- 
tem can be expressed in a form analogous to Eq. (5). 



^GDT(r) = £o(p)e'"^L'tAGDT(p, z), 



(7) 



and the optical potential defined as 

3 Try 



V(r): 



ki((x)i^-(Do + ir) 



n(r). 



(8) 



(9) 



In principle, the solution obtained from GDT pro- 
vides a good approximation only in the region where the 
optical potential is nonvanishing [19]. Nevertheless, it 
provides also a good approximation for calculating the 
scattering cross-sections using the formulae derived in 
Ref . [19] in the limit when a ^ \, Equation (8) can not 
be considered as a generally correct asymptotic form of 
the field. It is worth realizing, however, that it has a 
similar form to the solution of the Maxwell equation 
obtained with a slowly varying envelope approximation 
for the incident laser beam, with a polarization source 
given by a classical distribution of two-level atoms with 
a density n(r) [20]. In this sense, the optical potential 
can be interpreted as the polarization related to a two- 
level medium (with F = y). 

Both results as given above can be put in a more 
transparent form, namely the near field object plane 
field distribution can be written as 



£o(p) = ^o(p)2r(p). 



(10) 



where Sq (p) is the probe field amplitude, and assumed 
to be a Gaussian beam ^ exp(-pVpo^) with width 
Po = 1.5 mm, and ^(p) is the transmission function of 
the cloud. By doing so, we have expressed the scattering 
through the condensate as the transmission over an opti- 
cal element (BEC) whose transmittance is given by 
^(p) (we have neglected the free propagating phase 
factor e'^^^ here). This approximation requires the slowly 
varying scaled near field distribution to attain a stable 
asymptotic form at distance Za of the order of the size of 
the condensate, i.e., much smaller than the distance in 
which the image plane is located. Mathematically, this is 
equivalent to 



STip) ~ ijjip, z.) - tA(p, ^X 
a^ Za « observation distance scale (di, d2,f). 
In the on-shell approximation, we obtain 



2ros(p)-i + E ^^Bi(k,) 

1=0 ^ 



(11) 



r 



X sine deJo(kLpsme)Pi(cose) 



(12) 



with 
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where we have approximated cos^ - 1 ~ in the near 
field. In order to make specific calculations we need to 
know the atomic density profile n(r). In the following 
we shall use a spherically symmetric Gaussian profile 
given by 



From Eq. (13) it follows that the BEC is centered at 
r = 0. The near field patterns around the center as given 
by Eqs. (6) and (16) are presented in Figs. 1 and 2. 



n(r) = 



N 



(iTTSiy 



exp 



2a' 



(13) 



Such a profile is in fact quite appropriate for the case of 
atoms with a negative scattering length [9, 12]. For the 
case of a positive scattering length the density profile 
can be determined from the solutions of the (mean-field) 
Gross-Pitaevski equation (Ref. [21], see also Refs. [9, 
13]). As we have demonstrated in Ref. [13], however, 
using a Gaussian profile with an appropriately adjusted 
width leads to quantitatively similar results for scattered 
light. 

Using the density profile given by Eq. (13) we obtain 



(14) 



yl = 7Vy3V277e-^^^ lu.aikla^ilk^a) , 



where // denotes a modified Bessel function. Within the 
framework of the GDT the transmission function is 
given by 




= exp 



- 1 



Tgdt (p) 

7eff 



o}]^-(x>Q + tr 



exp - 



'2a' 



(15) 



(^> 



We can analogously obtain an expression for the scaled 
near field 



»Agdt(p, z) = exp 



_ ^"7eff expf— ^ 

0}]^- (x>Q + tr \ 2a^ 



^V\+Erf{z)] 



(16) 



where Erf denotes the error function. We stress once 
again that more general density profiles, including the 
asymmetric ones as well as those obtained from mean 
field Gross-Pitaevskii equation for BEC can be easily 
implemented into the present framework [13]. We will 
not discuss these issues here, since we expect such pro- 
files lead to qualitative and quantitative features of the 
scattered radiation very similar to those obtained with 
the Gaussian profile Eq. [(13)]. 



.</% 




Fig. 1. Near field distribution obtained from the on-shell approxima- 
tion given by Eq. (5) for 6 = 6y = {lir) 1.5 X 10^ MHz. The rest of 
the parameters is given in the text. The (maximum) on-axis (p = 0) 
optical density is about 0.5; (a) the amplitude, (b) the phase of 
(/^os(p, z) respectively. 
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Fig. 2. Same as in Fig. 1, except that obtained from the GDT given 
by Eq. (15). 



As we can see, already in a distance of the order of a 
few a beyond the condensate, stable asymptotic field 
patterns get established. Thus we can approximate the 
object plane to be also at r = 0. It should be pointed out 
that the near field patterns in Figs. 1 and 2 differ quite 
significantly. The GDT result exhibits a decrease of the 
amplitude and a phase shift close to the optical axis (for 
p — 0) as z grows. The result of the on-shell approxima- 
tion, on the other hand, varies in this regime much more 
slowly. In Fig. 3 we compare the results for the transmis- 
sion function as given by Eqs. (12) and (15). Here, the 
off-axial variation of the GDT result (dashed line) is less 



Fig. 3. The transmission functions obtained from the on-shell 
approximation Eq. (1 1) (sohd hne) and from the GDT Eq. (13) (dotted 
line). Both the phase and the amplitude are displayed for d = 67. 



significant than the variation of the result of the on-shell 
approximation (solid line). We attribute these differ- 
ences to the fact that the GDT describes in fact more 
accurately the gradual absorption and related refraction 
along the optical axis. The on-shell approximation 
reconstructs well the field inside or after the medium 
only, when the absorption has already taken place. As 
we shall see below, the differences between the two 
approaches decrease essentially in the process of optical 
imaging. 

In calculations we have used the following values of 
the parameters: A = 0.7 |uim (0.78 fxm for ^^Rb and 0.67 
imm for 'Li) , A^ = 2000, r=y= (2 77) 2.5 Mhz, and the 
trap size a = 2 |uim. The results enable us to make 
qualitative and to some extent quantitative comparisons 
with the experimental results obtained for ^'Rb and 'Li 
atoms respectively [1, 3]. 



3. Imaging Formation Process 

A (coherent) imaging system converts the near field 
pattern E^ip) in the object plane into a new distribution 
El (p) in the image plane. In this study the optical setup 
involves a set of compound (ideal) lens free of abbera- 
tions or any other kind of anomalies. The lens has a focal 
length /, and is located in a distance di away from the 
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object plane (i.e., from the atomic cloud), and in a dis- 
tance 6?2 away from the image plane (CCD camera). If 
the object is in the focus, then according to elementary 
optics 6?2 is determined by \ld\ + 21 d2 = l/f. The lens are 
accompanied by an aperture of diameter D in front. 

It is reasonable to assume that the whole imaging 
process is paraxial, and can be treated using the Fresnel 
diffraction theory [22] . The magnification of the above 
described system is M± = d2/d\. Since a magnifying 
optical system is not translationally invariant, the calcu- 
lation of Ei(p) from E^ip) has to proceed in two steps. 
First, we consider a magnification of the object neglect- 
ing the inversion of the image. Then we blur the image 
with an on-axis impulse response function (or, in an- 
other words, the Green's function) of the imaging sys- 
tem [22]. 

The derivation of the impulse response function is 
straight forward by following the propagation of the 
phase front of the field. An on-axis impulse in the object 
plane is propagated in free space in the Fresnel limit 
over a distance di to the face of the lens. It is then 
multiplied by a spatial filtering profile of the pupil func- 
tion due to the aperture, as well as an ideal phase retar- 
dation due the lens. This is followed again by a free 
space paraxial propagation up to the image plane. The 
impulse response function in such a situation is [22] 



^. .. ^ 2J,{7TD\p-pVXd2) 



(17) 



where 



Go = 7TD^/{4X^did2), and J^ is the Bessel function of 
order 1. The spatial resolution of this impulse response 
function in the image plane is ps ~ 1- llkdJD. 
Finally, the image can be constructed from 



^i(p) 



"Mi J 



dp'£o(pVMJGp(p-p') 



= ^ J dp'^o(pVMOSrfpyMOGp(p - p') . (18) 



where the pref actor \IM\ is from the reduction of image 
intensity due to magnification. 

The values of the parameters involved in the imaging 
process were chosen by us in such a way so that they 



correspond to the values reported in the Rice experi- 
ment [3] with / = 12 cm, D = 3 cm, d\ = 17 cm, 
and 6?2 = 40. 8 cm. The magnification of the system is 
therefore M±_ = 2.4, and the resolution of the impulse 
response function is p^ = 11.61 |ULm in the image 
plane, which corresponds to a resolving power of 
ps = ps/Mi = 3.7 |ULm in the object plane. From these 
numbers, we conclude that any structures of order less 
than ps in the object plane would correspond to a point 
source as far as the imaging properties are concerned. 



4. Results and Discussions 

We have used Eq. (1 8) to study the diffraction pattern 
formation. The numerical calculation can be done in 
two ways, either as a direct two-dimensional integration 
of Eq. (18) in terms of x' and j' [p' = (x\ j')], or by 
utilizing the Fourier transform technique. In the latter 
case, fast Fourier transform (FPH") algorithms can be 
used to compute the transforms for E^ip) and Gp(p) 
separately before the inverse of their product is made to 
yield the result. Both EFT and explicit integration meth- 
ods were tested by us. 

Using Babinet's principle [22], we can also calculate 
the shadow image E'l as generated by the complemen- 
tary transmission function ST'' = l-ST. This provides an 
alternative numerical scheme which is more suitable, 
since now the integration in Eq. (18) is basically limited 
to the region of the cloud. 

In Fig. 4 we present the field distribution in the image 
plane. As we mentioned, imaging reduces the differ- 
ences between the on-shell approximation and the GDT 
results. Both approaches give now qualitatively similar 
results. In particular, the phase of £'i(p) grows from 
negative values to zero for p of the order of p^, and 
remains practically constant then. The amplitude of the 
field exhibits characteristic oscillations. 

In actual CCD image measurements one measures 
the relative absorption coefficient, which can be defined 
as 



^(p) = 1 



i^i(p)r 



\E,(p) + ^ (pf 



(19) 



where E[(p) + Et(p) represent the image of the probe 
without the presence of the cloud. 

The relative absorption coefficients for the on-shell 
(solid line) and GDT (dashed line) approximations is 
shown in Fig. 5. This figure can be directly compared to 
Fig. 2 of Ref. [3]. The qualitative agreement in both 
cases is satisfactory. Quantitatively, the GDT result 
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5. Conclusions 

We have studied formation of diffraction patterns in 
the near-resonant imaging of trapped cold atomic sam- 
ples. We showed that spatial imaging similar to the 
spectral measurement of the coherent Rayleigh scatter- 
ing [13] can probe the detailed information of the 
trapped atomic clouds, and may provide direct evidence 
of the presence of the Bose-Einstein condensate with 
improved imaging optics. We have used the on-shell 
approximation and the Glauber's generalized diffraction 
theory to study light scattering off BEC and to construct 
near field distributions. Optical imaging theory was 
used to construct the field distributions in the image 
plane. The results compare very well with the recent 
experiments, and should be useful for the experimenters 
involved in BEC research. 

Acknowledgments 



Fig. 4. The image plane field distribution obtained from the on-shell 
approximation (solid line) and from the GDT (dotted line). Both the 
phase and the amplitude are displayed for 8 = 67. 
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Fig. 5. The relative absorption as given by Eq. (18) calculated from 
the on-shell approximation (solid line) and the GDT (dotted line). For 
8 = 67. 



is closer to experimental data, giving very similar values 
of the relative absorption coefficient. We attribute this 
result to the fact that GDT seems to describe absorbtion 
and refraction within the medium in a more accurate 
manner. The oscillation period in Fig. 5 for both results 
of the GDT and the on-shell approximation is, however, 
slightly shorter than that observed by Bradley et al. [3]. 



We thank Drs. E. Cornell, M. Holland and W Ketterle 
for many helpful discussions. The work of L. You is 
supported by the NSF through a grant for the Institute 
for Theoretical Atomic and Molecular Physics at 
Harvard University and Smithsonian Astrophysical 
Observatory. 



6. References 

[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, 
and E. A. Cornell, Science 269, 198 (1995). 

[2] K. B. Davis, M.-O. Mewes, M. R. Andrews, N. J. van Druten, D. 
S. Durfee, D. M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 
3969 (1995). 

[3] C. C. Bradley, C. A. Sackett, J. J. Tollett, and R. G. Hulet, Phys. 
Rev. Lett. 75, 1687 (1995). 

[4] For recent references see W. Petrich, M. H. Anderson, J. R. 
Ensher, and E. A. Cornell, Phys. Rev. Lett. 74, 3352 (9195); 
K. B. Davis, M. O. Mewes, M. A. Joffe, M. R. Andrews, and 
W. Ketterle, ibid, 5202 (1995); C. S. Adams, H. J. Lee, N. 
Davidson, M. Kasevich, and S. Chu, ibid, 3577 (1995), and 
references therein. 

[5] A. L. Fetter and J. D. Walecka, Quantum Theory of Many- 
Particle Systems, McGraw-Hill, New York (1971); Bose Einstein 
Condensation, A. Griffin, D. W. Snoke, and A. Stringari, eds., 
Cambridge University Press, Cambridge, (1995); S. T. Beliaev, 
Sov. Phys. JETP 7, 289 (1958) [Zh. Eksp. Teor. Fiz. 34, 417 
(1958). 

[6] J. R. Gardiner, R. A. Chne, J. D. Miller, D. J. Heinzen, 
H. M. J. M. Boesten, and B. J. Verhaar, Phys. Rev. Lett. 74, 3764 
(1995). 

[7] A. J. Moerdijk, W. C. Stwalley, R. G. Hulet, and B. J. Verhaar, 
Phys. Rev. Lett. 72, 40 (1994); R. Cote, A. Dalgarno, and M. J. 
Jamieson, Phys. Rev. A 50, 399 (1994); E. R. I. Abraham, 
W. I. McAlexander, C. A. Sackett, and R. G. Hulet, Phys. Rev. 
Lett. 74, 1315 (1995). 



581 



Volume 101, Number 4, July-August 1996 

Journal of Research of the National Institute of Standards and Technology 

[8] I. F. Silvera and J. T. M. Walraven, Spin -Polarized Hydrogen, in 

Progress in Low Temperature Physics, D. F. Brewer, ed., 

Elsevier, The Netherlands, (1986); J. T M. Walraven, in 

Proceedings of the Fourteenth International Conference on 

Atomic Physics, Boulder, Colorado, 1994, D. J. Wineland, 

C. E. Wieman, and S. J. Smith, eds., AIP Press, New York 

(1995); T J. Greytak and D. Kleppner, in New Trends in 

Atomic Physics, G. Grynberg and R. Stora, eds., North- 
Holland, Amsterdam (1984) p. 1125. 
[9] P A. Ruprecht, M. J. Holland, K. Burnett, and M. Edwards, 

Phys. Rev. A 51, 4704 (1995). 
[10] H. T C. Stoof, Phys. Rev. A 45, 8398 (1992). 
[11] M. Lewenstein and L. You, Phys. Rev. A 53, 909 (1996). 
[12] G. Baym and C. Pethick, Phys. Rev. Lett. 76, 6 (1996) 
[13] L. You, M. Lewenstein, R. J. Glauber, and J. Cooper, Phys. Rev. 

A 53, 329 (1996); this paper also contains an extensive list of 

references on light scattering. 
[14] B. Svistunov and G. Shlyapnikov, Zh. Eksp. Teor. Fiz. 97, 821 

(1990) [Sov. Phys. JETP 70, 460 (1990) ]; ibid. 98, 129 

(1990) [71,71 (1990)]. 
[15] H. D. Pohtzer, Phys. Rev. A 43, 6444 (1991). 
[16] J. Javanainen, Phys. Rev. Lett. 72, 2375 (1994). 
[17] (a) M. Lewenstein, L. You, and J. Cooper, in Proceedings of the 

Int. Conf. Quantum Optics HI — Szczyrk 1993, M. Kolwas and 

J. Mostowski, eds., special issue of Acta Phys. Polon. A 86, 173 

(1994) ; (b) L. You, M. Lewenstein, and J. Cooper, Phys. Rev. 

A 50, R3565 (1994). 
[18] J. R. Taylor, Scattering Theory: The Quantum Theory of Non- 

relativistic Collisions, Wiley, New York (1987). 
[19] R. J. Glauber, in Lectures in Theoretical Physics, W. E. Brittin 

and L. G. Dunham, eds.. Vol. I, Interscience, New York (1959) 

p. 315. 
[20] L. Allen and J. H. Eberly, Optical Resonance and Two-Level 

Atoms, Dover, New York (1987). 
[21] M. Edwards and K. Burnett, Phys. Rev. A 51, 1382 (1995). 
[22] B. E. A. Saleh and M. C. Teich, Chap. 4 in Fundamentals of 

Photonics, Wiley, New York (1991). 



About the authors: Li You is a National Science Foun- 
dation postdoctoral fellow at Harvard University. He 
conducts theoretical research at the Institute for Theo- 
retical Atomic and Molecular Physics (ITAMP). Maciej 
Lewenstein is a senior theorist at the Service de 
Photons, Atomes, et Molecules of the Commissariat a 
I 'Energie Atomique (Saclay, France). He conducts 
research in several domains of atomic, molecular, and 
optical physics, such as interactions of matter with 
intense laser fields, quantum optics, and the quantum 
field theory of atoms and photons. 



582 



